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ABSTRACT 


Following  success  in  other  areas  as  a  remote  sensor,  the  spacebome  microwave 
image  radars  are  assuming  a  notable  position  in  the  problem  of  ship  detection  for  civilian 
and  military  purposes.  This  work  will  discuss  the  strong  and  weak  points  of  Synthetic 
Aperture  Radar  (SAR)  when  used  for  ship  detection.  First,  the  thesis  gives  a  brief 
description  of  SAR  fundamentals,  image  processing  and  the  parameters  for  ship 
detection.  Second,  the  actual  techniques,  limitations,  errors  and  some  models  used  for 
ship  detection  are  described.  Finally,  using  a  well-known  and  reliable  ship  detection 
model  (Vachon  et  al.  1997),  tested  in  the  Canadian  Ocean  Monitoring  Workstation  and  in 
some  validation  field  programs,  the  new  generation  of  spacebome  SARs,  mainly 
RADARSAT  2,  are  analyzed  for  ship  detection  capabilities.  During  the  analysis 
parameters  like  wind  velocity,  wind  direction  related  to  the  antenna,  satellite  incident 
angle  and  Number  of  Looks  are  changed  to  study  their  influence  on  ship  detection. 
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I.  INTRODUCTION 


A.  MOTIVATION 

Having  approximately  7,480  km  (4650  miles)  of  coastline  and  200  nautical  miles 
of  territorial  sea,  the  total  marginal  sea  area  under  the  supervision  of  the  Brazilian  Navy 
totals  the  impressive  amount  of  approximately  2,773,000  kn?.  Maritime  control  and 
surveillance  of  this  huge  area  requires  all  types  of  tools.  To  fulfill  this  task  in  time,  with 
reliable  accuracy,  airborne  and  spaceborne  microwave  imaging  sensors  emerge  as  a 
valuable  instrument.  The  advantages  of  this  kind  of  tool  include  a  regular  repeating 
observation  cycle,  images  free  of  cloud  interference,  high  ground  resolution  (up  to  3m  in 
RADARSAT-2)  and  reliability.  Synthetic  Aperture  Radars  (SAR)  have  these  essential 
features  and  offer  a  good  trade-off  between  cost  and  benefit.  In  the  future,  use  of  SAR 
will  be  supported  by  a  commercial  remote  sensing  imaging  boom.  The  new  commercial 
SAR  systems  can  supply  images  for  specific  requests  with  competitive  cost  and  a 
potential  for  integration  within  modular  ship  monitoring  systems.  An  increasing  number 
of  missions  are  being  planned  and  executed  as  outlined  in  Table  1. 

This  thesis  focuses  specifically  on  three  missions  that  will  fly  the  next  generation 
SAR  sensors  designed  to  fulfill  various  applications  including  ship  detection.  These 
sensors  will  be  launched  between  2001  and  2003  and  include  the  following:  the  Japanese 
Space  Agency  Advanced  Land  Observing  Satellite  with  the  Phase  Array  type  L-band 
Synthetic  Aperture  Radar  (ALOS/PALSAR);  the  European  Space  Agency  Environmental 
Satellite  with  an  Advanced  Synthetic  Aperture  Radar  (ENVISAT/ASAR)  and  the 
Canadian  Space  Agency  Radar  Satellite  2  (RADARS AT  2). 
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Table  1. 


Planned  SAR  satellites  missions  [Updated  from  Olsen  et  al.  (2000)] 


Satellite 

Description 

1  Key  radar  parameters 

(Owner/Operator) 

Swath  width 

Incidenc 

spatial 

Modes 

e  angles 

res. 

ENVISAT  (USA) 

ENVISAT  will  be  the  next 

56-105  Km 

15-45** 

30  m 

Image 

Launch 

satellite  to  be  launched  by  the 

56-105  km 

15-45° 

30  m 

ALTPOL 

jan  2002 

European  Space  Agency.  The 

405  km 

17-42° 

150  m 

Wide 

C-band  ASAR  instrument  will 
be  the  first  dual-polarised 
orbital  SAR  with  a  long  term 
presence  in  space  [3]. 

405  km 

17-42° 

1000  m 

Global 

RADARSAT-2 

This  mission  is  designed  to  be 

20-50° 

28  m 

Standard 

(MDA/CSA) 

yet  another  step  closer  to  a 

20-45° 

28  m 

Wide 

Launch  2003 

fully  commercial  mission. 

10-20° 

28  m 

Low  Inc 

Most  of  the  investment  is 

50-60° 

20  m 

High  Inc 

carried  by  the  Orbital 

37-48° 

9m 

Eine 

Sciences  Subsidiary,  Mac 

500  km 

20-50° 

100  m 

ScanSAR  Wide 

Donald  Dettwiler,  with  the 

300  km 

20-46° 

50  m 

ScanSARNarrow 

Canadian  Space  Agency  as  an 

20-41° 

Std  QuadPod 

anchor  customer.  The  satellite 

25  km 

28  m 

Pine  QuadPod 

will  carry  a  C-band  SAR  that 

25  km 

30-41° 

9m 

Triple  Pine 

supports  viewing  to  either  side 

50  km 

30-50° 

3  m 

Ultra  Pine  Wide 

of  the  spacecraft  track,  high 

20  km 

30-40° 

3  m 

Ultra  Pine  Narrow 

and  low  resolution  modes  and 
polarimetric  capability. 

10  km 

30-40° 

ALOS  (NASDA) 

Primarily  a  land  observing 

300  km 

18-48° 

10-20  m 

Image/Dual  pol 

Launch  2003 

mission,  operating  at  L-band. 

modes 

The  wide  swath  mode 

100  m 

ScanSAR 

provides  an  imaging  geometry 
that  will  provide  good  target 
to  clutter  ratios. 

Terras  AR 

Terras  AR  is  a  commercial 

X-band: 

(DLR/DASA/MMS/BN 

mission  planned  by  a 

10  km 

20-55° 

1.6  m 

Spotlight 

SC)  [4] 

consortium  of  British  and 

20  km 

25-45° 

3  m 

Stripmap 

German  industry.  Current 

100  km 

20-45° 

15  m 

ScanSAR 

plans  indicated  a  two-satellite 

L-band: 

system  with  X  and  L-band 

40  km 

20-35° 

9m 

capability. 

60  km 

20-45° 

9  m 

Quadpol 

200  km 

20-45° 

30  m 

Stripmap 

ScanSAR 

COSMO/SKYMED 

This  is  presently  a  mission 

3  km 

<45° 

1  m 

Spotlight 

concept  developed  under 

30  km 

<45° 

3  m 

Stripmap 

funding  from  the  Italian  Space 
Agency.  The  SAR  system 
considered  is  the  so-called 
SAR2000,  an  X-band  radar 
capable  of  stripmap,  spotlight 
and  ScanSAR  imaging. 

120  km 

<45° 

25  m 

ScanSAR 

LightSAR 

This  mission  has  been  planned  for  some  time  by  NASA,  where  an  industry  partner 

(NASA/Industry) 

has  been  sought.  Currently,  discussions  are  being  held  with  oil  industry,  to  determine  how  a 
mission  can  be  designed  to  meet  that  industry  sector’s  requirements  as  well  as  NASA’s  science 

1  goals.  The  proposed  radar  is  an  L-band  system  with  moderate  resolution. 

Earth  Watch 

1  The  ESA  Earth  Watch  programme  is 

intended  to  be  ESA’s  contribution  towards  a  I 

(ESA/Industry) 

European  Earth  observation  and  monitoring  system,  in  cooperation  with  private  industry  and 

organizations  such  as  EUMETSAT.  Industry  proposals  for  mission 

concepts  have  been 

solicited,  and  several  SAR  missions  are  included.  Eurther  selections  have  not  been  confirmed 

1  to  date. 
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This  thesis  presents  a  technical  review  of  one  of  these  spaceborne  SARs  -  the 
Canadian  RADARSAR  2.  This  choice  is  due  to  the  availability  of  commercial  imagery, 
access  to  public  information  on  the  sensor,  the  success  of  the  previous  mission 
(RADARSAT-1),  the  set  of  technical  improvements,  and  the  suitability  for  the  Brazilian 
Navy  Project.  Issues  like  ship  tracking  and  target  recognition  that  result  from  ship 
detectability  improvements  will  be  mentioned  but  not  pursued  in  detail  in  this  thesis. 
Expected  results  with  new  polarization  improvements  will  also  be  mentioned  but  not 
analyzed  in  detail.  This  thesis  will  describe  SAR  operations  and  the  physics  of  the  ship 
detection  problem. 

B.  THESIS  OUTLINE 

The  first  chapter  describes  the  motivation  and  focus  of  this  thesis.  Chapter  II 
presents  some  necessary  SAR  fundamentals  and  background  needed  to  understand  the 
ship  detection  problem.  Constraints,  limitations  and  explanation  of  system  modes 
precedes  the  SAR  imaging  process  in  order  to  present  a  complete  view  of  data  acquisition 
through  final  image  creation.  Chapter  III  deals  q)ecifically  with  all  facets  of  the  ship 
detection  problems,  limitations,  errors  and  techniques.  The  fourth  chapter  is  the  thesis 
core.  All  the  results  of  quantitative  improvements  for  the  next  family  of  spaceborne 
SARs  are  included  in  Chapter  IV.  Finally,  Chapter  V  summarizes  the  results  and 
conclusions  focusing  on  SAR  applicability  to  ship  detection. 
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n.  BACKGROUND 


A.  THE  SYNTHETIC  APERTURE  RADAR  (SAR)  FUNDAMENTALS 

1.  Basics 

SAR  systems . 

A  Synthetic  Aperture  Radar  (SAR)  is  a  Pulse- Doppler  imaging  radar  that 
transmits  short  radio  pulses  in  the  microwave  region  of  the  spectrum.  A  SAR  obtains 
resolution  in  range  through  travel  time  measurement,  and  uses  Doppler  shift  of  the 
backscattered  signal  to  obtain  resolution  in  azimuth  as  shown  in  Figure  1.  The  azimuth 
resolution  is  obtained  using  the  concept  of  the  synthetic  antenna  or  synthetic  aperture 
(Figure  2).  A  very  long  antenna  is  synthesized  by  moving  a  small  antenna  along  a 
convenient  path  on  a  platform  flight  path  or  azimuth  direction.  The  returning  Doppler- 
shifted  backscattered  signal  (amplitude  and  phase)  is  digitally  processed  to  obtain 
information  in  the  azimuth  direction.  The  combination  of  time-resolved  range  and 
frequency-resolved  azimuth  results  in  a  SAR  image  of  target  backscatter. 

The  SAR  systems  create  a  large  amount  of  data  necessitating  extensive  processing 
to  produce  images  with  the  required  resolution.  This  process  is  usually  performed  and 
stored  in  ground  stations.  In  the  early  days  of  SAR  observations,  the  processing  was  done 
optically,  but  due  to  the  huge  increase  in  computational  power,  all  processing  is  currently 
done  digitally. 

The  SAR  technique  was  extended  to  use  two  antennas  leading  to  the  generation  of 
three-dimensional  images  of  the  illuminated  surface  (Graham,  1974).  This  technique  is 
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usually  referred  to  as  across  track  Interferometric  Synthetic  Aperture  Radar  (IFAR); 


however,  IFAR  is  not  the  subject  of  this  thesis. 


Figure  1.  Azimuth  and  Range  direction. 

[From  http://southport.ipl.nasa.gov/desc/imagingradarv3.html1 


■  S'j^ilhe^c  Aperture  Length,  L " 


Figure  2.  Azimuthal  Synthetic  Aperture 

[From  http://southport.ipl.nasa.gov/desc/imagingradarv3.html 
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2. 


SAR  theory 


The  foundation  of  SAR  theory  is  the  synthesized  antenna  aperture  to  increase  the 
along  track  or  azimuthal  resolution.  The  SAR  theory  is  described  in  Curlander  and 
McDonough  (1992),  and  Elachi  (1982).  Figure  3  in  three  dimensions  and  Figure  4  in  two 
dimensions  summarize  all  terminology  and  geometry  for  SAR.  A  full,  detailed 
description  of  SAR  theory  is  beyond  the  scope  of  this  work  but  some  basics  are  needed  to 
better  understand  the  ship  detection  problem  described  below. 
a)  Azimuthal  resolution 

The  well-known  result  shown  below  is  described  by  Curlander  and 
McDonough  (1992)  among  others  and  began  with  the  azimuth  resolution  for  Side- 
Looking  Aperture  Radar  (STAR).  The  resolution  is  independent  of  the  range  and  the 
radar  wavelength,  and  improves  with  a  smaller  aperture,  given  by 

rap=PA  ~  [1] 

where  da  is  the  antenna  real  aperture  in  azimuthal  direction. 

The  corresponding  resolution  expression  for  an  orbital  SAR  returns  the  complete 

result 


Pa 


Re  d, 

Re  +  h  2 


[2] 


where  Re  is  the  earth  radius  and  h  is  the  platform  altitude  (Oliver,  1998). 
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b)  Range  resolution 

The  minimum  separation  between  two  resolvable  points  in  SAR  can  be 


defined  as  the  minimum  range  separation  of  two  points  that  can  be  distinguished  by  the 
system  as  separate.  For  SAR  sensors 

AR,  cx 

=  =— 

sinr]  2sinr| 

or  after  pulse  compression  techniques 


where  c  is  the  light  speed,  AR^  is  the  resolution  of  the  slant  range,  r\  is  the  incident  angle 
between  the  radar  beam  and  normal  to  the  earth’s  surface  and  x  ^  is  pulse  duration. 

c)  Design  considerations  (constraints  and  limitations) 

Below  is  a  summary  of  the  highlights  of  the  critical  design  aspects  of  a 
SAR  as  an  emitting  and  receiving  sensor.  Other  technical  construction  constraints  related 
to  SARs  are  beyond  the  scope  of  this  work. 

(1)  Choice  of  Pulse  Repetition  Frequency  (PRF)  or  Pulse 

constraints. 


The  correct  pulse  duration  is  needed  to  achieve  a  reasonable  range 

resolution  and  a  sufficient  echo  signal  to  noise  ratio  (SNR).  Chirp  pulses  and 

compression  techniques  are  used  to  achieve  the  optimum  PRF.  This  is  a  severe  constraint 

in  designing  practical  SARs  due  to  the  compromise  between  improving  azimuth 

resolutbn  and  providing  a  wider  swath  width.  Therefore,  ignoring  the  earth’s  effect  and 

orbit  curvature  for  simplicity,  we  have  from  Robinson,  (1985)  the  following: 
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c 


[5] 


where  PRF  is  Pulse  Repetition  Frequency,  W  is  the  swath  width,  V  is  the  platform 
velocity,  ?  is  the  angle  between  the  slant  direction  and  vertical  (incident  angle).  The 
minimum  PRF  limitation  is  easier  to  visualize  in  terms  of  antenna  theory  because  in  those 
terms  what  determines  the  minimum  PRF  is  the  distance  between  successive  array 
elements  in  the  real  antenna.  If  this  distance  is  greater  the  X/2,  the  so-called  grating  lobes 
(replications  of  the  arrays’  main  lobe  at  increasing  intervals  on  either  side  of  the  main 
lobe)  will  be  produced. 

(2)  Sidelobes  Minimization. 

Stimson,  1998  states. 

Performance  of  a  synthetic  array  radar  may  be  degraded  by  both  range  sidelobes 
due  to  pulse  compression  and  the  sidelobes  of  the  synthetic  array.  The  sidelobes 
affect  the  radar  maps  in  two  different  ways.  First,  the  peaks  of  the  stronger 
sidelobes  may  cause  a  string  of  progressively  weaker  false  targets  to  appear  on 
either  side  of  a  strong  target.  Second,  the  combined  power  of  all  sidelobes- called 
integrated  sidelobe  return- together  with  noise  tends  to  fog  or  wash  out  the  detail 
of  the  maps. 

(3)  Motion  Compensation 

Platform  velocity  and  orbit  are  not  perfectly  constant  as  SAR 
theory  assumes.  The  whole  SAR  concept  is  based  on  very  slight  signal  phase  changes 
over  comparatively  long  periods.  Therefore,  it  is  essential  that  any  motion  during 
short  periods  be  accounted  for.  A  complete  study  of  this  error  source  is  beyond  the 
goal  of  this  thesis,  but  the  reader  is  referred  to  AGARD-LS-182  (1992). 
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(4) 


Antenna  size  limitations 


The  reduction  in  antenna  size  is  limited  by  the  operational 
necessity  to  generate  a  sufficiently  powerful  signal  in  order  to  obtain  enough  returned 
signal.  A  minimum  value  of  the  vertical  antenna  dimension  is  needed  to  focus  the 
beam  into  the  desired  swath  width.  These  two  design  constraints  are  other  practical 
limitations  to  SAR  system  resolution  improvements. 

3.  SAR  operating  modes 

There  are  three  operational  modes  of  SAR  systems:  STRIPMAP,  SCAN  AND 
SPOTLIGHT,  pictorially  sketched  in  Figure  5  (Franceschetti,  1999).  These  STRIPMAP 
is  the  most  popular  which  was  developed  and  adapted  into  the  following  modes: 

•  Squint  forward  and  backward  -  an  array  that,  although  it  loses  some 
azimuthal  resolution,  can  execute  applications  like  targeting/detection 
ahead  of  and  behind  the  platform,  and  backscattering  analysis  of  the 
properties  of  the  illuminated  surfaces  with  respect  to  the  azimuth  angle. 

•  Doppler  beam  sharpening  provides  high  quality,  continuously  updating 
maps  of  large  expansions  of  the  surface  with  the  penalty  of  reducing  cross¬ 
range  resolution  distance  Multilook  mapping.  This  mode  improves  map 
quality  by  averaging  out  the  scintillation  of  the  radar  returns. 

•  The  SCAN  mode  tremendously  increases  the  range  swath  dimension  at  the 
expense  of  azimuthal  resolution  while  stepping  the  antenna  beam  to 
neighboring  sub- swaths  in  range  direction. 
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Figure  3. 
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SAR  terminology  and  geometry  1  [From  AGARD-LS- 182, 1992] 
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Figure  4.  SAR  terminology  and  geometry  2  [From  AGARD-LS-182,  1992] 
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•  The  SPOTLIGHT  mode  is  steered  during  overall  acquisition  time  to 
illuminate  the  same  area.  A  huge  gain  in  azimuthal  resolution  is  traded  for 
a  loss  in  coverage  (small  illuminated  area  along  the  sensor  path). 


Figure  5.  SAR  operational  modes  [From  Franchescetti,  1999] 
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B. 


SAR  IMAGE  PROCESSING 


Using  Stewart,(1985),  the  SAR  image  processing  produces  high  spatial  resolutbn 
in  the  azimuthal  direction  and  high  spatial  resolution  in  range.  Processing  the  radar  signal 
in  essence  consists  of  taking  a  scattering  time  history  at  each  range  and  convolving  it 
with  the  point-target  response  to  produce  a  line  in  the  image.  Using  a  continuous 
approach,  Franceschetti  (1999)  states  the  two-dimensional  processed  SAR  image 
expression  as  follows: 


r')  =  |jY(x.r)sin  c 


—  ix'-x) 
Ax 


sm  c 


dxdr 


[6] 


A 

where  y{x,r)  represents  a  two  dimensional  reflectivity  scene  pattern  including  the  phase 
factor.  Symbols  x  and  r  are  the  azimuthal  and  range  direction  respectively. 

1.  Geometric  Resolutions 

From  Franceschetti  (1999),  also. 


Simply  speaking,  geometric  resolution  is  the  ability  of  the  system  to  localize 
nearby  objects.  More  precisely,  the  resolution  length  is  the  minimum  spacing 
between  two  objects  that  are  detected  as  separate  entities,  and  are  therefore 
resolved. 


Both  results  for  range  and  azimuthal  were  already  presented  in  [1]  and  [3]. 
a)  Focused  x  Unfocused 

An  “unfocused”  SAR  processing  system  is  one  which  attains  its  along- 
track  resolution  by  simple  frequency  filtering  of  the  Doppler  waveform.  Higher  azimuthal 
resolutions  are  only  possible  in  a  focused  SAR  system,  which  can  account  for  the 
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variation  in  Doppler  frequency  of  a  target  as  it  passes  through  the  footprint  resulting  in 
the  along-track  resolution  limit  of  half  the  antenna  length. 

2.  Radiometric  Resolutions 

Franceschetti,  (1999)  defines  radiometric  resolution  as  “  measure  of  the  ability  of 
the  system  to  discriminate,  or  resolve,  areas  of  different  scattering  properties.  These  are 
described  by  the  reflectivity  pattern  (y(x,r))  of  the  illuminated  surface.” 


Franceschetti,  (1999)  also  states  that  changes  in  the  reflective  pattern  are  related 
to  two  different  processes,  one  macroscopic  (image  intensity  modulation  proportional  to 
Y(x,r)^  due  to  surface  shape  and/or  its  electromagnetic  parameters  as  a  function  of  space 
coordinates)  and  the  other  microscopic  (essentially  phase  change  due  to  surface 
roughness).  A  conventional  definition  of  radiometric  resolution  given  by  Brooks  and 
Miller,  (1979)  and  Franceschetti,  (1999)  covers  both  microscopic  and  macroscopic 
process  follows: 


^rd  =10*  log 


V 


[7] 


where  pis  the  mean  value  and  o  is  the  standard  deviation  of  the  image  intensity 

referring  to  a  distributed  target  of  constant  reflectivity.  This  states  that  adjoined  areas 
with  different  intensities  can  be  resolved,  provided  that  their  difference  is  larger  than 
Ard. 


3.  Radiometric  and  Geometric  Distortions 

In  practice,  many  factors  can  affect  image  quality  by  causing  signal  amplitude  or 
phase  modulation.  The  former  is  a  consequence  of  antenna  failure  in  stabilization 
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creating  a  change  in  the  antenna  point  angle.  Phase  modulation  is  a  consequence  of  an 
unpredicted  change  in  the  range  between  the  object  and  the  antenna,  giving  rise  to  an 
unknown  receiving  signal,  which  affects  azimuthal  pulse  compression. 

Oliver,  (1998)  summarizes  that  the  amplitude  modulation  leads  to  radiometric 
distortions  and  increased  side  lobes,  and  phase  modulation  leads  to  image  defocus, 
geometrical  distortions  (foreshortening,  layover  and  shadow),  radiometric  distortions  and 
increased  sidelobe  levels.  Image  defocus,  radiometric  distortion,  and  increased  sidelobe 
levels  restrict  target  detection  and  recognition  capability.  Geometrical  distortions  prevent 
direct  pixel-to-pixel  comparison,  which  is  relevant  to  multichannel  registration,  target 
change  detection,  and  direct  image  comparison.  Motion  Compensation  (MOCO),  made 
by  accurate  accelerometers  and  inertial  navigational  units  mounted  on  the  antenna,  can 
reduce  amplitude  modulation  effects.  However,  after  this  processes  some  residual 
unknown  phase  modulation  still  exists. 

4.  SAR  Signal  Statistics 

In  the  case  of  a  SAR  scene,  the  only  way  to  describe  its  roughness  and  density  of 
scatters,  for  surface  and  volumetric  scattering  respectively,  is  in  terms  of  statistical 
parameters.  Data  characterization  throughout  a  probability  density  function  (PDF)  is  a 
matter  of  compromise  between  simple  models  with  few  degrees  of  freedom  and  complex 
models  with  many  degrees  of  freedom.  Some  well-known  PDF  with  two-degrees  of 
freedom  like  Weibull,  lognormal,  and  the  K  distribution  are  used  to  statistically  represent 
a  SAR  image.  We  will  focus  on  the  K-distribution  since  it  is  the  PDF  used  in  our  ship 
detection  model. 
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5.  Speckle 

The  fluctuation  in  the  receiving  signal,  known  as  fading,  caused  by  the  diverse 
nature  of  the  return  echo,  generates  a  noise  like  grainy  appearance  on  SAR  imagery 
called  speckle.  It  influences  the  ability  to  estimate  image  properties  and  thus  is  central  to 
information  retrieval  from  individual  SAR  images.  Speckle  has  a  noise  like  appearance, 
but  it  is  not  noise.  An  electromagnetic  measurement  can  be  useful,  for  example,  in  SAR 
interferometry.  It  can  also  be  understood  as  an  interference  phenomenon  in  which  the 
principal  source  of  the  noise  like  quality  of  the  observed  data  is  the  distribution  of  the 
phase  terms. 

C.  RELEVANT  TOPICS  IN  SHIP  DETECTION 

1.  Radar  Cross  Section  (RCS) 

The  RCS  describes  the  backscattering  property  of  the  target  and  depends  on  its 
size,  shape  and  orientation  as  well  as  on  wavelength  and  polarization  of  the  incident 
signal  (Skolnik,  2001).  From  the  definition  of  RCS  (o)  the  cross  section  of  an  object  that 
scatters  equally  in  all  directions  (isotropic)  is  equal  to  its  projected  area. 

a)  Normalized  Radar  Cross  Section  (NRCS) 

The  quantity  o°.  Radar  Cross  Section  per  Unit  Area,  provides  a 
normalized  parameter  that  can  be  used  to  describe  RCS  of  the  Sea  independent  of  the 
radar  illumination. 

The  normalized  radar  cross  section  (NRCS)  is  given  by 

a°  =  a/A  [7] 

where  A  is  the  area  of  a  smooth  surface  that  corresponds  to  the  mean  land  or  sea  surface 

area  contained  within  the  radar’s  cell  resolution.  The  illuminated  area  is  a  complicated 
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function  of  radar  system  parameters  and  is  generally  predictable  as  is  described  below  in 
Section  3. 

2.  Ship  Radar  Cross  Section 

Ship  reflectance  is  a  complicated  problem  without  a  general  solution.  Skolnik 
(2001)  proposed  an  empirical  relation  between  displacement  in  tons  and  ship  radar  cross 
section.  Vachon  et  al.  (1997)  using  the  Skolnik  relationship  developed  the  following 
linear  regression  between  ship  weight  and  length  for  the  Bedford  Institute  of 
Oceanography  Fleet  and  for  some  ships  participating  in  the  MARCOT’95: 

a  =  D  =  O.OS.L^^^  [8] 

where  L  is  ship  length  and  D  is  displacement. 

Olsen  et  al.  (2000)  improved  the  above  formula  including  the  effect  of  incident 
angle  using  the  same  data  from  MARCOT’95 

a  =  D  =  R(0).O.O8.L^^^  [9] 

where 

R(e)  =  0.78+0.11.0  [10] 

All  these  formulations  show  that  ship  RCS  is  founded  onobserved  values  and  empirical 
methods. 

3.  Ocean  Radar  Cross  Section 

a)  Nature  of  o°  for  Ocean  and  Incident  Angle  Dependency 

The  NRCS  is  sensitive  to  the  character  of  the  sea  surface  and  is  highly 

variable.  It  can  change  as  much  as  10  dB  in  a  one- minute  interval  due  to  variations 

among  measuring  instruments  and  their  calibration  Long  (2001).  Therefore,  measurement 
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errors  and  uncertainties  tend  to  obscure  the  weak  functional  relationship  between  the 
wavelength  and  o°  for  most  incident  angles.  This  relationship  is  a  complicated 
function  of  incidence  angle.  For  instance,  at  low  incidence  angle  near  the  vertical,  a°  is 
highly  dependent  on  sea  state  and,  for  high  incident  angles,  a°  is  a  strong  function  of  A, . 
For  horizontal  polarization,  the  curve  of  a°  versus  the  grazing  angle  ^  6  for  low 
frequencies  is  expected  to  be  the  same  for  microwaves  except  for  a  larger  critical  angle. 
Therefore,  the  curve  for  horizontal  polarization  might  appear  as  in  Figure  6. 

b)  Dependence  of  <5°  on  Polarization 

Polarization  affects  o°  for  different  surface  conditions  or  for  different 
wavelengths.  From  Long  (2001)  makes  the  following  statements  that  hold  for  the  ratio 
between  a  vv  (NRCS  for  vertical  transmitted  and  vertical  received  polarization)  and  a  hh 
(NRCS  for  horizontal  transmitted  and  horizontal  received  polarization): 

1 .  The  ratio  increases  with  an  increase  in  wavelength, 

2.  The  ratio  decreases  with  an  increase  in  sea  roughness,  and 

3.  c°hh  can  exceed  o°vv  for  heavy  seas  and  small  depression  angles. 

It  should  not  be  erroneously  assumed  from  the  above  that  under  very  rough 
conditions  a°vv  and  a°hh  are  about  equal  for  all  incident  angles. 

c)  Dependence  of  o°  on  Wind  and  Sea 

Some  of  the  sea  surface  characteristics  known  to  influence  various 
features  of  radar  echos  include  the  period  and  shape  of  the  waves,  the  wave  height,  wind 
ripples,  and  the  presence  or  absence  of  whitecaps  and  spray.  However,  observations 

1  0°  incidence  angle  from  nadir  (vertical)  is  equal  to  90°  grazing  angle  measured  from  horizontal  plane 
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indicate  that  average  and  medium  values  of  c°  increase  as  the  sea  becomes  rougher  for 
all  grazing  angles  except  for  those  near  vertical  incidence.  Vachon  et  al.  (1997)  uses  an 
empirical  polarization  ratio  constructed  by  Unal  et  al.  (1991)  to  transform  from  modeled 
values  of  a°vv  to  values  of  a°hh  necessary  for  the  ship  length  computation. 


Figure  6.  Shape  of  o°  versus  Grazing  Angle  curves  for  horizontal  polarization 
[From  Long,  2001] 
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m.  SHIP  DETECTION  PROBLEM 


A.  SHIP  DETECTION  APPROACHES 

1.  Ship- Generate  Surface  Waves  (Wake  detection) 

Ship  wake  patterns  exhibit  various  characteristics  in  SAR  images  under  diverse 
environmental  conditions.  Wake  structures  in  an  image  can  be  classified  in  three 
categories.  First,  surface  waves  generated  by  the  ship;  second,  turbulent  wakes  or  vortices 
containing  relatively  persistent  but  non-propagating  currents  visible  in  SAR  imagery 
through  the  interaction  of  ambient  waves  with  these  surface  currents  (Lynden,  1988). 
Third,  internal  waves  generated  by  the  ship  that  are  visible  due  to  interaction  with  the 
short  surface  waves. 

Ship  wakes  that  are  observed  in  SAR  imagery  are  created  by  two  mechanisms. 
First,  variations  of  short  waves  (centimeters- scale)  appearing  as  bright  narrow  “V”  wakes 
are  directly  observed  due  to  Bragg  scattering  in  images;  second,  longer  waves  (decameter 
scale)  forming  the  classical  Kelvin- wake  system  are  observed  as  they  modulate  the 
formation  of  Bragg  scattering  waves.  Standard  ship  wake  detection  techniques  exploit 
the  first  subcategory  due  to  better  pixel  intensity  (bright)  response  in  SAR  images.  Some 
techniques  (Rey  et  al.  (1990),  Copeland  et  al.  (1995)  and  others)  have  been  developed  to 
use  the  Radon  transform^  in  extracting  linear  features  that  are  correlated  with  ship  wakes. 

Many  specialists  consider  point  target  approaches  a  more  efficient  “pattern 
recognition”  technique  than  wake  detection.  Wake  detection  techniques  are  very  useful, 
after  a  candidate  vessel  is  chosen,  in  order  to  confirm  or  reject  the  target  as  ship. 

2  A  commonly  used  technique  for  finding  lines  in  an  image.  Very  robust  in  the  presence  of  high  levels 
of  noise  but  lacking  in  differencing  between  short  and  long  lines  or  end  points  of  short  lines. 
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Furthermore,  in  key  applications  like  fisheries  monitoring  the  vessels  may  often  travel 
too  slowly  to  generate  significant  wake  signature  (Olsen,  2000).  This  work  will  present  a 
performance  analysis  of  new  spacebome  SARs  using  a  point  target  detection  approach. 

2.  Point  Target  Detection  (Ship  detection) 

The  main  task  of  point  target  detection  is  to  detect  a  single  or  a  small  cluster  of 
pixels  with  the  higher  backscatter  intensity  from  a  ship  within  the  background  scene 
intensity.  Generally,  the  background  ocean  has  a  lower  backscatter  intensity  than  ships. 
The  ocean  backscatter  varies  significantly  with  wind  speed  and  wind  direction  relative  to 
the  SAR  view  angle,  SAR  beam  incidence  angle,  and  polarization  of  the  returned  power. 
The  challenge  of  ship  detection  is  unambiguously  identifying  the  small  number  of  bright 
ship  pixels  among  the  variable  intensity  of  the  background  population  of  pixels. 

3.  Ship  Detection  Techniques 

An  intuitive  approach  to  ship  detection  is  to  divide  the  scene  into  small  frames 
assuming  a  single  ocean  backscatter  value  within  each  frame  that  determines  a  detection 
threshold  for  candidate  vessels.  The  simplest  method  is  to  set  the  threshold  intensity  at  n 
standard  deviations  above  the  image  mean  intensity.  One  variation  of  this  technique  uses 
two  windows  inside  the  image  with  different  sizes.  The  windows  are  moved  across  the 
entire  image  using  the  larger  window  to  compute  the  ocean  background  statistics  while 
the  smaller  window  is  used  for  the  actual  target  search. 

Some  more  sophisticated  parametric  approaches  fit  a  PDF  (usually  as  a  K- 
distribution  )  to  the  data  and  estimate  the  threshold  based  on  the  uncertainty  level  desired. 
One  non-parametric  method  uses  a  probabilistic  Neural  Network  approach.  All  these 
methods  are  designed  to  reach  a  balance  between  detection  rate  and  the  occurrence  of 
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false  alarms.  To  achieve  this,  some  analysis  methods  are  applied  such  as  contrast 
measures,  homogeneity  tests,  morphological  filters,  wake  analyses,  and  combinations  of 
these. 

a)  The  Constant  False  Alarm  ( CFAR )  method 

The  CFAR  method  sets  a  threshold  for  radars  operating  in  a  sea-clutter 
environment  based  on  a  PDF  for  the  image  ocean  background  intensities,  integrated  until 
a  required  significant  level  is  reached.  This  significant  level,  Tjc,  corresponds  to  a  CFAR 
of  l-Tjc  for  the  intensity  PDF 

p{x)dx  [8] 

Two  main  problems  arise  from  this  method.  The  first  problem  is  to  find  a  PDF  that  fits  an 
image  intensity  histogram.  Second,  F  must  be  computed.  Vachon  (1997)  uses  r|c  = 
0.995  and  p(x)  is  defined  as  a  K-distribution  for  various  order  parameters. 
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IV.  NEXT  GENERATION  SHIP  DETECTION  EXPECTED 

RESULTS 


The  main  question  to  be  answered  by  a  Ship  Detection  Model  is  what  size  of 
vessels  can  be  detected  in  SAR  imagery  and  how  reliable  is  the  detection  under  various 
observational  conditions  (Olsen,  2000).  Secondarily,  the  model  can  indicate  whether 
information  like  type,  speed  and  heading  can  be  retrieved  from  SAR  imagery. 

The  key  operational  requirements  for  the  next  generation  spaceborne  SARs  that 
are  needed  to  improve  ship  detection  are: 

•  Increase  in  resolution  to  be  applied  in  target  detection; 

•  increase  in  swath  coverage  to  allow  widespread  ship  monitoring  and  surveillance; 

•  larger  incident  angles  (>  35°)  to  reduce  ocean  clutter  and  increase  ship 
backscattering; 

•  greater  coverage  frequency  in  combination  with  higher  resolution  and 
programmable  modes  for  selected  areas. 

Almost  all  of  these  desirable  aspects  will  be  flying  on  the  next  spaceborne  SARs.  The 
goal  of  this  thesis  is  to  apply,  in  a  reliable  model,  the  new  specifications,  mainly  for 
resolution  and  incident  angle,  and  retrieve  the  general  performance  behavior  for  different 
satellites  and  their  different  modes. 
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A. 


RELEVANT  NEXT  GENERATION  SATELLITE  PARAMETERS 


This  section  gives  a  brief  description  of  all  three  satellite  and  sensor 
characteristics  and  the  launch  schedule  over  the  next  2-3  year  window.  Details  for 
RADARSAT2,  ENVISAT/ERS2  and  ALOS/PALSAR  SAR  systems  are  summarized  in 
the  following  sections 


RADARSAT-  2 


1. 


Table  2  presents  RADARSAT-2  values  for  swath  width,  incident  angles.  Number  of 
Looks,  azimuthal  and  range  spatial  resolutions,  multipolarization  capability  and  others 
characteristics,  like  orbit  parameters  and  coverage  access  using  500km  swath  width. 
Figure  7  shows  the  coverage  pattern  of  all  operational  modes  and  Table  3  shows 
RADARSAT-2  inovations. 


Figure  7.  RADARSAR  2  Modes  [From  CCRS/RADARSAT-  2  site] 
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Table  2.  RADARSAR  2  characteristics  [From  CCRS/RADARSAT-  2  site] 


Standard  Beam 
Modes 

Nominal 

Swath 

Width 

Incidence 

Angles 

Number 
of  Looks 

Approx. 

Resolution 

Wide 

100km 

20-50 

1x4 

25m  X  28m 

Low  Incidence 

1x4 

25m  X  28m 

High  Incidence 

1x4 

40m  X  28m 

Line 

70km 

50-60 

1x4 

20m  X  28m 

ScanSAR  Wide 

50km 

37-48 

1x1 

10m  X  9m 

ScanSAR  Narrow 

500km 

20-50 

4x2 

100m  X 

100m 

Standard  Quad  Pol. 

300km 

20-46 

2x2 

50m  X  50m 

Line  Quad  Pol. 

25km 

20-41 

1x4 

25m  X  28m 

Triple  Fine 

1 

11m  X  9m 

Ultra- fine  Wide 

3x1 

1  Im  X  9m 

Ultra- fine  Narrow 

20km 

30-40 

1 

3m  X  3m 

ORBIT  CHARACTERISTICS 


Altitude  (average) 

798  km 

Inclination 

98.6° 

Period 

100.7  minutes 

Ascending  Node 

18:00  hrs 

Sun-  synchronous 

14  orbits  per  day 

Repeat  Cycle 

24  days 

COVERAGE  ACCESS  USING  500km  SWATH  WIDTH 


North  of  70°  H 

Daily 

North  of  48°  H 

Every  1-2  days 

Equator 

Every  2-3  days 
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Table  3.  RADARSAT-2  innovations  [From  Sanden  (2001)] 


INNOVATIONS 

BENEFITS 

GPS  receivers  onboard 

+60- meter  real-time  position 

information 

10  ms  delay  between  imaging 
modes 

Faster  mode  changes 

Yaw- steering  for  zero- Doppler 
shift  at  beam  center 

Facilitates  image  processing 

Higher  downlink  power  density 

3 -metre  minimum  size  antenna  on 
ground  allowing  station  portability 
and  Lower  "cost  of  entry"  for  new 
ground  stations 

3 -meter  ultra- fine  resolution 

Highest- resolution  commercially 

available  SAR 

Left- and  right- looking  capability 

Faster  revisit  time,  2000  km 
accessibility  swath  and  Routine 
Antarctic  mapping  available 

Fully  polarimetric  imaging  modes 

Enhanced  capabilities  for  various 
applications 

Solid-state  recorders  for  onboard 
image  storage 

Higher  reliability,  faster  image  access 
and  Simultaneous  record  and 
downhnk 

2.  ENVISAT/ASAR 

ASAR  has  five  mutually  exelusive  modes  of  operation  whieh  could  be  classified 
in  two  categories:  Global  mission  subdivided  in  to  Global  monitoring  and  Wave  mode 
that  have  a  low  data  rate  and  hence  an  operation  capability  up  to  100%  of  the  orbit,  and 
regional  mission  subdivided  in  narrow  swath  modes,  image  mode  and  alternating 
polarization  mode.  Table  4  summarizes  some  operational  parameters  and  Table  5  the 
mode’s  characteristics. 
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Table  4. 


ENVISAT/ASAR  Operational  Parameters  [From  ENVISAT  Internet  site] 


Instrument 

Parameters 

Image  Mode 

Alternating 

Polarization 

Wide  Swat 
h 

Global 

Monitoring 

Wave  Mode 

Swath  width 

Up  to  1 00  km 

Up  to  100 
km 

>  400  km 

>  400  km 

5  km  vignette 

Operation 

time 

Up  to  30  min  per  orbit 

Rest  of  orbit 

Data  Rate 

Up  to  100  Mbit/s 

0.9  Mbit/s 

Power 

1365  W 

1395  W 

1200  W 

713  W 

647  W 

Table  5.  ENVISAT/ASAR  Modes  [From  ENVISAT  Internet  site] 


ASAR 

Swathes 

Range 

Incidence 

Angle 

Swath  Far 

Width  Range 
[km] 

Near  Incidence 
Angle 

IS1 

108.4  -  109.0 

14.1  -  14.4 

22.2  -  22.3 

IS2 

107.1  -  107.7 

18.4  -  18.7 

26.1  -26.2 

IS3 

83.9  -  84.3 

25.6  -  25.9 

31.1  -31.3 

IS4 

90.1  -90.6 

30.6  -  30.9 

36.1  -36.2 

IS5 

65.7  -  66.0 

35.5  -  35.8 

39.2  -  39.4 

IS6 

72.3  -  72.7 

38.8-39.1 

42.6  -  42.8 

IS7 

57.8  -  58.0 

42.2  -  42.6 

45.1  -45.3 

3.  ALOS/PALSAR 

The  Advanced  Land  Observing  Satellite  with  a  Phased  Array  type  L-band 
Synthetic  Aperture  Radar  (PALSAR)  is  the  Japanese  next  generation  spaceborne  SAR. 
PALSAR  uses  Lband  frequencies,  different  from  RADARSAT  and  ENVISAT/ASAR, 
and  has  land  observation  as  a  major  task  and  therefore  has  a  variable  beam  elevation.  The 
development  of  the  PALSAR  is  a  joint  project  between  NASDA  and  Japan  Resources 
Observation  System  Organization  (JAROS).  Its  major  specifications  are  listed  in  Table  6. 
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Table  6.  ALOS/PALSAR  specifications 


[From  ALOS/PALSAR  site] 


Observation  Mode 

Eine  Resolution  Mode 

ScanSAR  Mode 

Frequency 

L-band 

Polarization 

HH  or  VV  (option:  HV  or  VH) 

Spatial  Resolution 

10m  (21ooks)/20m  (4  looks) 

100m 

Swath  Width 

70km 

250-360km(3  -  5scans) 

Off- nadir  Angle 

18  -  48deg. 

S/A 

25dB 

NE?* 

-25dB 

B.  SHIP  DETECTION  MODEL  APPLIED  ON  RADARSAT-1 

Vachon  et  al.  (1997)  develops  a  model  to  study  RADARSAT-  1  (C  band  HH 
polarization)  performance  for  ship  detection  using  the  various  beam  modes.  The  model 
was  validated  with  data  acquired  during  the  RADARSAT-1  ship  detection/validation 
field  program.  Minimum  detectable  ship  length  is  used  as  a  Figure  of  Merit  (FOM)  for 
ship  detection  using  RADARSAT-1  standard  modes  and  ERS-1  modes.  Validation  of  the 
model’s  key  assumptions  are  provided,  including  those  concerned  with  the  hybrid  C-band 
HH  polarization  ocean  radar  cross  section  model,  the  image  probability  density  function 
and  the  ship  radar  signature  model.  This  model  is  being  updated  by  the  Canadian  Center 
of  Remote  Sensing  (CCRS)  and  was  applied  and  tested  in  the  Ocean  Monitoring 
Workstation  -  a  system  developed  by  Satlantic,  Inc.  with  technical  and  financial  support 
of  CCRS,  The  Department  of  Fisheries  and  Oceans  (DFO),  Canadian  Coast  Guard 
(CCG),  the  Department  of  National  Defense  (DND)  and  the  Canadian  Space  Agency 
(CSA).  The  OMW  is  designed  to  provide  operational  users  of  marine  data  with  near-real 
time,  value  added  ocean  information  derived  from  RADARSAT-1  SAR  images  (Vachon, 
2000). 
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1.  Model  Assumptions 

The  model  is  based  on  three  assumptions: 

(1)  the  SAR  image  probability  density  function  for  ocean  scenes  can  be 
described  by  the  compound  K-distribution  model  for  sea  clutter; 

(2)  the  radar  cross  section  of  the  ocean  can  be  estimated  given  wind  speed  and 
beam  geometry  of  C-Band  radars; 

(3)  ship  size  is  directly  related  to  the  radar  cross  section  of  the  ship. 

The  model  was  developed  using  ERS-1  SAR  data  prior  to  RADARSAT-1  launch  and 
was  subsequently  validated  with  RADARS  AT- 1  data. 

2.  SAR  Ocean  Scene  Image  Statistics 

The  compound  K-distribution  model  is  an  empirical  model  for  sea  clutter  that  is 
applicable  to  high-resolution  radars  that  resolve  fine  structure  on  the  sea  surface. 
Gaussian  distributions  do  not  fit  for  this  kind  of  processes  and  the  K-distribution  does  so 
due  to  the  following  observations  from  Rey,  1996: 

The  envelope  time  history  of  a  demodulated  signal  from  an  individual  range  cell 
shows  fast  fluctuations  from  pulse  to  pulse  that  is  modulated  by  an  underlying 
structure. 

The  fast  fluctuation  component,  or  ‘speckle’,  decorrelates  from  pulse  to  pulse 
with  frequency  agile  signals  but  is  correlated  for  between  5  and  10  ms  with  fixed 
frequency  signals.  This  implies  a  large  number  of  scatters  within  each  illuminated 
patch.  By  applying  the  central  limit  theorem,  the  ‘speckle’  can  be  approximated 
by  the  Rayleigh  distribution. 

The  Slower  modulation  component  is  the  local  sea  clutter  mean  level,  which  has  a 
long  temporal  decorrelation  period  and  is  not  affected  by  frequency  agility.  It 
characterizes  the  mean  level  variation  of  clutter,  including  clutter  ‘spikes’,  and  its 
intensity  is  found  to  be  gamma  distributed; 

The  gamma  distribution  is  found  in  this  model  to  be  the  first  term  in  a  Laguerre 
series  expansion  of  the  modulation  PDF.  The  implication  is  that  the  non- Gaussian 
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nature  of  sea  clutter  seems  to  be  due  to  more  to  scatter  bunching  by  the  sea  wave 
structure  than  anything  else. 


The  proposed  distribution  to  describe  the  SAR  PDF  of  image  intensity  is  the 
multi- look  intensity  K- distribution  (Oliver,  1993): 


p{l)=j 


f  Lt)/  ^ 

(L+f))/2 

9 

1/2 

[<’>) 

r(L)r(t)) 

(</>j 

[9] 


where  p(I)  is  the  probability  of  the  image  intensity  I,  </>  is  the  mean  image  intensity,  L 
is  the  number  of  statistically  independent  looks,  v  is  an  order  parameter  for  the  intensity 
modulation  process,  F  is  the  gamma- function,  Kv-l  is  the  modified  Bessel- function  of  v-L 
order. 


Using  the  above  PDF,  we  can  now  compute  the  critical  intensity  using  [8]  for  any 
significant  level,  r|c.  To  retrieve  the  same  results  as  Vachon  et  al.  (1991)  we  use  0.995 
significance  level,  a  trapezoidal  numerical  routine  for  integration  (0.005  and  0.01 
intervals),  and  image  mean  </>  =  1. 

3.  RADARSAT- 1  SAR  parameters 

The  RADARSAT- 1  parameters  are  summarized  on  Table  7.  The  values  are  for 
nominal  (near  beam  center)  parameters  for  the  21  modes  of  the  RADARSAT- 1  SAR. 
ERS- 1  values  are  also  given  for  comparison. 

4.  Radar  Cross  Section  of  the  Ocean 

A  physically-based  model  for  the  a°  relevant  to  RADARSAT- 1  C-Band  HH 
polarization  is  not  available.  Vachon  et  al.  (1997),  using  Gower  et  al.  (1993),  based  on 
Snoeij  et  al.  (1991),  incorporates  empirical  results  from  Unal  et  al.  (1991)  for  backscatter 
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polarization  ratio  (PR“).  Table  8  lists  the  conversion  factors  to  convert  from  a  VV  to  HH 
Ocean  Radar  Cross  Section  polarization  model. 

Table  7.  RADARSAR-1  parameters  [From  Vachon  (1997)] 


Beam  Mode 

Swath  [km] 

Looks  (L) 

0[deg] 

Pazim] 

Prdm] 

o°„,[dB] 

Dchirp  [km] 

SI 

100 

3.1 

23.5 

27 

24.2 

-23 

15.8 

S2 

100 

3.1 

27.5 

27 

20.4 

-21 

13.6 

S3 

100 

3.1 

33.5 

27 

25.3 

-23 

11.4 

S4 

100 

3.1 

37.0 

27 

23.4 

-23 

10.5 

55 

100 

3.1 

39.0 

27 

22.1 

-22 

10.0 

S6 

100 

3.1 

43.5 

27 

20.3 

-24 

9.1 

57 

100 

3.1 

47.0 

27 

19.1 

-23 

8.6 

R 

50 

1 

55.5 

8.3 

-23 

10.1 

E2 

50 

1 

40.5 

8.4 

7.9 

-22 

9.7 

E3 

50 

1 

42.5 

7.6 

-23 

9.3 

E4 

50 

1 

44.5 

8.4 

7.3 

-24 

9.0 

F5 

50 

1 

46.0 

8.4 

7.1 

-23 

8.6 

W1 

150 

3.1 

25.5 

27 

33.8 

-22 

14.6 

W2 

150 

3.1 

35.0 

27 

24.6 

-22 

11.0 

W3 

150 

3.1 

42.0 

27 

20.8 

-25 

9.4 

EHI 

75 

3.1 

50.5 

27 

18 

-25 

8.1 

EH2 

75 

3.1 

51.5 

27 

17.7 

-25 

8.0 

EH3 

75 

3.1 

53.5 

27 

17.3 

-25 

7.8 

EH4 

75 

3.1 

55.5 

27 

16.8 

-23 

7.6 

EH5 

75 

3.1 

57.0 

27 

16.6 

-23 

7.5 

EH6 

75 

3.1 

55.5 

27 

16.4 

-23 

7.3 

ELI 

170 

3.1 

16.5 

27 

39.1 

-22 

22.2 

ERS-I 

100 

4.9 

23.0 

30.7 

38.4 

-24 

14.2 

This  table  includes  the  following  parameters  and  modes:  Sstandard;  Ffine;  W-  wide;  EH- 
extended  high  incidence;  EL -extended  low  incidence;  L  is  the  number  of  independent  looks;  0  is  the 
incidence  angle  paz  is  the  azimuth  resolution;  prg  is  the  ground-range  resolution;  o^ne  is  the  noise- 
equivalent  clutter  level;  Dchirp  is  the  ground-projected  chirp  length. 


Table  8.  Vertical  to  Horizontal  C-Band  (5.3  GHz)  Polarization  Ratios  [From  Unal 

et  al.  (1991)] 


Inc.ang/ 

Wind(m/s) 

20° 

o 

o 

45 

2 

1.05 

2.07 

4.78 

4 

0.79 

2.37 

5.25 

6 

0.65 

2.57 

5.47 

8 

0.56 

2.70 

5.59 

10 

0.51 

2.80 

5.66 

12 

0.49 

2.88 

5.71 

14 

0.49 

2.95 

5.75 
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The  well-developed  C-band  VV  (vertical/vertical)  polarization  model  (o°vv)  for 
ERS-1  scatterometer  wind  retrieval  (the  C  band  Models-CMOD  series),  was  used  here  in 
conjunction  with  Table  8  ratios  to  produce  the  HH  (horizontal/horizontal)  RCS  of  the 
Ocean  (a°hh)  as  showed  below 

a°hh  =  a°vv- PR“  [10] 

5.  Radar  Cross  Section  of  Ships 

Based  on  Skolnik  (2001),  Vachon  et  al.  did  a  linear  regression  between  ship 
weight  and  length  for  the  Bedford  Institute  of  Oceanography  (BIO)  fleet  and  some  ships 
participating  in  the  military  exercise  MARCOT’95,  resulting  in  the  following  empirical 
relation: 

Cship=  D  =  [11] 

where  Oship  is  the  ship  radar  cross  section,  ‘D’  is  the  ship  displacement  and  ‘L’  ship 
length. 

The  model  can  then  determine  the  ship  detectability  by  a  SAR  based  on  a  ship 
scale  parameter.  Knowing  the  critical  intensity  (T)  of  the  relevant  PDF  from  [9]  and  the 
ocean  normalized  cross  section  (a°)  from  [10],  the  model  can  determine  the 
corresponding  minimum  point  target  RCS  for  detection  at  a  chosen  probability  level  as 

[12] 

where  paz  is  the  satellite  azimuthal  resolution  cell  size  and  prg  is  the  ground  range 
resolution  cell  size. 
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6.  Ship  Length  computation 

Combining  [11]  and  [12]  we  have  the  minimum  detectable  ship  length  as  a 
function  of  the  Ocean  RCS,  resolutions  and  critical  intensity,  as  follows: 

L  =  (I,o“p„P„/0.08)’"  [13] 

Applying  this  results  to  RADARSAT-1  data,  Vachon  et  al.  (1997)  produced  results  for 
SI,  S3  and  S7  RADARSAT-1  modes  and  ERS  to  test  the  model,  considering  -23  dB  as 
the  RADARSAT-1  noise  equivalent  normalized  radar  cross  section  (7ne)  with  clutter 
noise  ratio  (CNR)  =  a°/aNE  and  One  =  -24  dB  for  ERS. 

7.  Vachon  et  al.  Ship  Detection  Model  Results 

Some  results  for  this  model  for  RADARSAR  1  mode  SI,  S3,  S7  and  ERS  SAR 
are  presented  in  Figure  8.  Each  plot  shows  the  minimum  detectable  ship  length  as  a 
function  of  wind  speed  using  equation  [13]  where  solid  lines  represent  the  Upwind 
direction  (ie,  (|)=0,  the  largest  a°,  hence  the  worst  case  -  least  detectable  due  to  a  bright 
background)  while  dashed  lines  represent  the  Crosswind  direction  (ie,(l)=90°,  the  smallest 
a°,  hence  the  best  case).  The  three  curves  represent  variation  of  the  order  parameters  for 
K- distribution  from  top  to  bottom:  v=  4,  10,  8  (ie,  from  worst  to  best  case).  To  facilitate 
comparison  among  the  various  RADARSAT-1  beams  and  modes,  Vachon  et  al.  (1997) 
defines  a  FOM  for  ship  detection  using  Upwind  (j)  =  0°  direction,  wind  velocity  U  =  10 
m/s,  and  order  parameter  v  =  4.  The  resulting  plot  is  shown  in  Figure  9. 
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8. 


Vachon  et  al.  results  validation 


A  validation  field  program  was  held  in  March/April  1996  to  provide  quantitative 
information  to  validate  the  ship  detection  model.  Some  of  the  elements  used  to  validate 
the  model  are  showed  in  Vachon  et  al.  (1997)  and  summarized  here.  They  include 
acquisition  of  wind  and  wave  data  from  two  buoys,  acquisition  of  nine  RADARSAT-1 
SAR  scenes  in  various  modes,  placement  of  ship  near  the  buoys  at  the  time  of 
RADARSAT- 1  SAR  passes,  and  surveillance  flights  over  the  area  to  identify  additional 
ship  targets  on  an  opportunity  basis. 

After  RADARSAT-1  calibration  and  SAR  data  processing,  Vachon  et  al.  (1997) 
concludes. 

We  validated  some  of  the  ship  detection  model  assumptions  using  our  calibrated 
RADARSAR  SAR  data.  In  particular,  we  show  that  there  is  excellent  agreement 
between  the  observed  RADARSAT- 1  SAR  cross  section  and  those  predicated  by 
the  hybrid  model.  Furthermore,  the  K- distribution  in  shown  to  be  a  suitable  PDF 
for  RADARSAT-1  SAR  ocean  images.  We  also  show  that  our  very  simple  ship 
length- dependent  cross  section  model  is  within  the  correct  order  of  magnitude 

C  SHIP  DETECTION  MODEL  APPLIED  TO  NEXT  GENERATION 

SPACEBORNE  SAR  PARAMETERS 

1.  Some  Considerations 

From  Vachon  et  al.  (1997),  we  have  some  important  statements  about  the  radar 
cross  section  of  the  ship.  They  are:  “...For  the  ship  considered,  we  find  that  our  simple 
relationship  is  within  20dB,  but  that  the  anomaly  tends  to  increase  with  increasing 
incidence  angle”,  and  more  “  ...Although  we  have  not  considered  ship  orientation  at  all 
in  this  analysis,  use  of  the  best  fit  line  to  our  data  allows  recovery  of  the  vessel  length  to 
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Figure  9.  Figure  of  Merit  from  Vachon  et  al.  (1997) 

within  20%  based  on  the  measured  ship  cross  section  alone.  It  is  apparent  that  our  simple 
ship  model  will  tend  to  underestimate  the  ship  detectability.” 

Due  to  the  good  results  presented  and  supported  by  Vachon  et  al.  (1997),  the 
model  is  used  here  to  describe  the  physical  processes  that  limit  ship  detection  and  to 
quantify  the  ship  detection  improvements  that  can  be  expected  in  the  next  generation 
spacebome  SARs. 

Unal’s  et  al.  (1991)  table  is  used  as  reference  and  including  only  values  linearly 
interpolated  inside  the  table  wind  range  (except  for  mode  S7  where  an  linear 
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extrapolation  to  47°  of  incident  angle)  The  computations  begin  at  2  m/s  wind  velocity  for 
to  two  reasons.  First,  due  to  a  lack  of  Unal  et  al.  (1991),  values  inside  this  interval  to 
confirm  our  values.  Secondly,  because  our  computed  curves  have  a  better  agreement  with 
Vachon  et  al.  (1997),  only  results  between  the  2  to  14  wind  speed  interval  are  included 
here. 


2.  Reproducing  Results  with  Vachon  et  al.  model: 

The  Vachon  et  al.  (1997)  model  is  reproduced  here  using  RADARSAT-2 
specifications  and  the  same  image  parameters  used  for  RADARSAT-1  except  for 
resolution.  All  considerations  taken  in  IV-C-1  were  generated  during  this  computation 
due  to  a  small  difference  in  length  values  and  a  more  linear  result  for  all  RADARSAT-2 
and  ERS  modes.  Figures  10,12,14  and  16  from  below  with  Figure  8. 

3.  Expected  Result  for  all  RADARSAT  1  modes  kept  in 
RADARSAT-2 

All  four  cases  resulting  from  [13]  show  the  same  ship  length  curve  behavior  as 
wind  speed  increases  for  the  same  fixed  values  used  by  Vachon  et  al.  (Figure  8).  The 
unchanged  beam  parameters  were  Number  of  Looks  (3.1),  each  mode  incident  angle,  and 
azimuthal  and  range  resolution.  Figures  11,13,15  and  17  show  the  cumulative  K- 
distributions  for  each  case  that  were  used  to  produce  an  critical  intensity  value  for  each 
mode. 

Table  9  shows  the  comparison  between  Vachon  et  al.  (1997)  results  and  the 
results  produced  for  this  thesis.  The  differences  in  minimum  detectable  ship  fength  are 
typically  less  than  5  m.  Therefore  the  reproduction  of  the  model  is  acceptable  for  the 
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main  purpose  of  this  thesis,  the  performance  prediction  for  next  generation  ship  detection 


spacebome  SAR. 


Figure  10.  SI  mode  ‘Wind  versus  Ship’  length  thesis  results 
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Figure  11.  SI  mode  Cumulative  PDF  with  Critical  Intensity  Values 


Figure  12.  S3  mode  ‘Wind  versus  Ship’  length  thesis  results 
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Figure  13.  S3  mode  Cumulative  PDF  with  Critical  Intensity  Values 
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Figure  14.  S7  mode  ‘Wind  versus  Ship’  length  thesis  results 
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Figure  15.  S7  mode  Cumulative  PDF  with  Critical  Intensity  Values 


Table  9.  Vachon  versus  thesis  length  differences  for  SI,  S3,  S7  and  ERS  modes 


Obs.:  Values  rounded  to  meter  and  values  computed  for  wind  =  2  m/s 


Modes 

Slv 

SI, 

Slv-Sl, 

S3^ 

S3, 

S3y-  S3t 

Worst 

27 

29 

+2 

13 

12 

-1 

Best 

20 

24 

-4 

10 

9 

-1 

Modes 

S7„ 

S7, 

S7y-  S7t 

ERSy 

ERSt 

ERSy.  ERSt 

Worst 

9 

4 

-5 

40 

38 

-2 

Best 

7 

3 

-4 

28 

30 

+2 

Worst  is  Upwind  with  Order  Parameter=4;  Best  is  Crosswind  with  Order  Parameter  =  inf 
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Figure  16.  ERS  ‘Wind  versus  Ship’  length  thesis  results 
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Figure  17.  ERS  mode  Cumulative  PDE  with  Critical  Intensity  Values 
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4.  Expected  Results  for  Fine  Mode  Resolution 

Figures  18,  19  and  20  show  the  results  using  the  RADARSAT-2  fine  mode 
resolutions.  The  mean  ineident  angle  for  each  fine  mode  and  the  released  Number  of 
Looks,  listed  in  Table  2,  were  used  for  each  mode.  All  the  Figures  show  the  expected 
minimum  ship  length  behavior  as  wind  speed  increases.  Values  with  crosswind  (best) 
and  upwind  (worst)  conditions  are  shown  with  three  different  K- distribution  order 
parameters.  Despite  the  disadvantage  of  lack  of  coverage  due  to  a  small  swath  width,  all 
expected  detected  ship  length  values  associated  with  new  high  resolutions  for  fine  modes 
(1.5  to  10  m)  will  lead  not  only  to  ship  detection  but  also  ship  recognition. 


Figure  18.  Ship  length  for  RADARS AR  2  Triple  Fine  mode 
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Figure  19.  Ship  length  for  RADARSAR  2  Ultra- Fine  Wide  mode 
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Figure  20.  Ship  length  for  RADARSAR  2  Ultra- Fine  Narrow  mode 
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5.  Ship  Length  versus  Resolution  Qualitative  Expected  Behavior 

Figure  21  shows  how  detectable  ship  length  changes  with  a  fixed  wind  velocity 
using  one  RADARSAR-2  mode  (Standard),  varying  with  resolution  for 
upwind/crosswind  (worst/best  case),  and  with  two  Number  of  Looks  values.  Three  major 
differences  are  noticeable.  First,  the  high  incidence  angle  detection  can  be  increased  in 
wide  swath  modes  as  much  as  five  times  compared  to  low  incidence  angle  observations 
(ScanSAR  narrow).  Next,  intermediate  resolution  modes  (Standard),  have  less  incident 
angle  influence  and  almost  the  same  influence  of  sea  RCS  and  wind  incident  angle. 
Finally,  at  high  resolution  all  three  variables  have  about  the  same  influence  on  final  ship 
length  computation. 

6.  Ship  Length  versus  Incident  Angle  Qualitative  Expected  Behavior 

Figure  22  shows  the  variation  of  detectable  ship  length  versus  incidence  angle  for 
a  fixed  wind  velocity  and  RADARSAR  mode  (standard  in  this  case)  varying  with 
resolution,  upwind/crosswind  (worst/best  case)  and  Number  of  Looks.  With  an  increasing 
incident  angle  the  ship  detectability  increases  by  a  factor  of  almost  2.5  for  this  specific 
case.  High  incident  angle  we  have  the  best  results  for  ship  detection. 

Figure  21  and  22  confirm  the  major  role  of  resolution  and  incident  angle  in  the 
Vachon  et  al.  (1997)  ship  detection  model.  The  effect  of  resolution  decreases  as  the 
incident  angle  increases.  Number  of  Looks  and  sensor- relative  wind  direction  have  a 
secondary,  but  not  disregardable  influence. 
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Figure  21.  Ship  Length  versus  Resolution  Results 


3S 

Incident  Anglej^^ 


ShHU  Laiiain  v&  &3#ellWe  Irwldant  Antfa  fat  0=1C  nVs,  Upwind  A  Cpo&&w^n(lfcw  ftADARSAT  1  $lendand  mode 


'  Upwli'HUC:rQss'^d  dlilfei 
^  Number  !>T  LMkfi  dinaienc 


Phi^.Rrt=9a  ^  =1 

—  Phi=0,Rea*99A  Loolta  *3.1 
+  Phi®O.Rea"7W  i  Lcoks  -1 
+  Phi=CI  .R«=7eW  A  L<?fl|is=3.1 
'  PhiH^  A  Looka  "1 

■  Phi=90'.RM=99‘  £,  Looks  =3.1 
PHi=30.Rcj^i=7aa  ^  Lcote  =1 
^  PhjjSO-Res-^rdd  a*  Looks  «3.1 


^  ResoAitlofi  ijffer-wicesm 


I  *  * 

a  60|-»  * 

-t  +  .4. 


+  + 


+  * 


+  +  +  +-, 


Ship  Length  versus  Incident  Angle  Results 
48 


Figure  22. 


7.  Practical  Application 

Given  the  wind  condition  as  in  Figure  23  at  ON,  34W  and  4N,  34W  (North 
Brazilian  Coast)  from  PIRATA  project,  we  can  use  the  Vachon  et  al.  (1997)  model  to 
estimate  ship  detection  performance.  For  instance,  using  the  same  parameters  and 
statistics  as  in  the  tested  model  for  the  SI  mode  case  on  RADARS  AT  2,  we  can  expect 
that  SAR  images  taken  under  the  same  conditions  as  in  Figure  10  (mode  wind  velocity  of 
6  m/s),  ship  lengths  from  32  to  43  (best  and  worst  case)  are  expected  to  be  detected.  This 
short  example  testifies  to  the  model’s  usefulness  in  analyzing  real  situations  setting  the 
best  parameters  for  ship  detection  with  RADARSAT  2. 


Figure  22.  Wind  histogram  from  PIRATA  site 
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V.  CONCLUSIONS 


The  ship  detection  model  outlined  by  Vachon  et  al.  (1997)  was  used  to  verify  and 
quantify  ship  detection  improvements  expected  in  the  next  generation  spacebome  SAR 
imagery.  This  thesis  focuses  on  the  influence  of  incidence  angle,  wind  speed  and 
resolution,  on  the  detectability  of  ships  in  SAR  imagery.  The  results  are  compared  to 
RADARSAT  2  parameters  to  test  the  improvement  expected  from  new  capabilities.  The 
impact  of  new  polarized  modes  in  ship  detection  is  not  tested  in  this  thesis. 

Vachon  et  al.  (1997)  uses  a  model  of  the  ocean  Radar  Cross  Section  (RCS)  as  a 
function  of  wind  speed  and  incidence  angle.  The  minimum  ship  length  that  can  be 
detected  against  the  background  ocean  radar  cross  section  is  determined  by  the  critical 
intensity  level  obtained  from  a  statistical  relationship  between  ship  size  and  RCS. 

The  minimum  detectable  ship  length  is  strongly  dependent  on  wind  speed  due  to 
the  changes  imposed  by  the  wind  on  RCS  of  the  ocean  and  on  the  critical  intensity.  For 
instance  in  RADARSAT-2  SI  mode,  with  3.1  Number  of  Looks,  with  23.5°  of  incident 
angle  with  99.5  %  of  confidence,  we  expect  winds  to  change  from  2  to  10  m/s  while  the 
detectable  ship  lengths  change  from  29m  to  54m  in  the  worst  case  and  from  24m  to  37m 
in  the  best  case.  For  areas  surveyed  by  the  PIRATA  project  on  the  North  Brazilian  Coast 
(ON,  34W  and  4N,  34W)  with  mode  wind  velocity  of  6  m/s,  minimum  detectable  ship 
lengths  from  32  m  to  43  m  (best  and  worst  cases)  are  expected  to  be  detected. 

The  minimum  detectable  ship  length  changes  dramatically  with  beam  incident 
angle.  The  increase  in  incidence  angle  reduces  Ocean  RCS  due  to  a  reduction  of  Bragg 
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scatter  to  the  sensor  and  the  increase  of  the  ship  RCS  due  to  an  increase  in  double  bounce 
scatter.  For  instance,  for  RADARS  AT  2  standard  mode,  with  10  m/s  wind  speed  in  the 
worst  case  (upwind  related  to  the  antenna),  increasing  the  incident  angle  from  20°  to  45° 
changes  the  detectable  ship  length  from  99m  to  20m  -  nearly  a  factor  of  five  change. 

The  minimum  detectable  ship  length  is  strongly  dependent  on  radar  resolution 
since  reduction  of  the  pixel  size  increases  the  probability  that  photons  from  that  pixel  are 
reflected  by  the  ship.  For  instance,  changing  modes  in  RADARS  AT  2  from  ScanSAR  to 
UltrafineNarrow,  thereby  increasing  resolution  by  a  factor  of  10^  in  area,  with  10  m/s 
wind  speed,  20°  of  incident  angle  and  Number  of  Looks  equal  1,  the  detectable  ship 
length  will  decrease  from  280m  to  15m  (worst  case)  -  almost  a  factor  20  change. 

The  best  ship  detection  results  are  obtained  with  RADARSAT  2  Fine  Resolution 
Modes.  The  minimum  detectable  ship  lengths  of  1.5  to  3.5  meters  for  fine  mode 
resolution  not  only  illustrate  the  improvement  of  ship  detection,  but  also  make  ship 
recognition  possible  given  the  right  conditions.  Additionally,  ship  recognition  will  be 
improved  by  cross  polarization. 

The  general  level  of  improvement  for  RADARSAT  2,  ALOS/PALSAR  and 
ENVI/ASAR  arises  from  the  new  resolution  modes  and  cross  polarization  Some  specific 
design  factors  should  be  mandatory  for  future  SARs  designed  specifically  for  ship 
detection.  High  incident  angles  orbits,  shorter  repeat  cycle,  multipolarization  and  large 
swath  width  at  high  resolution  are  among  the  main  SAR  characteristics  needed  for  future 
ship  detection  satellites. 
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APPENDICES 


A.  UNAL’S  TABLE  INTERPOLATION  CODE 

function  [Pr_Interpol_in_dB ]  =  Unals Int (V, theta, IntModeType ) 

o 

%  Interpolation  to  find  Unal ’ s  et  al .  tabulated  polarization  ratio 

o 

%  inputs 

%  V  -  the  wind  speed  in  m/ s 

%  theta  -  the  incidence  angle  in  degrees 

%  IntModeType  -  Matlab  interpolation  type,  ie,  'linear cubic' , etc 

%  output 

%  Pr_Interpol_in_dB  =  Polarization  Ration  Interpolated  in  decibels 

Beam_Inc_Tb=  [20  30  45  47]  ’ 

Spd_Wind_Tb=  [2  4  6  8  10  12  14] 

%  Unal ’ s  et  al . ,  table  values  in  dB  extrapoleted  until  47  degrees 

PRu=[  1.05  2.07  4.78  4.48 
0.79  2.37  5.25  5.25 
0.65  2.57  5.47  5.47 
0.56  2.70  5.59  5.59 
0.51  2.80  5.66  5.66 
0.49  2.88  5.71  5.71 
0.49  2.95  5.75  5.75]’; 

%Polari zat ion  Ratio  in  dB 

Pr_Interpol_in_dB  =  interp2  ( Spd_Wind_Tb,  Beam_Inc_Tb,  PRu,  V,  theta, ... 

IntModeType) ; 
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B.  CMOD-IFR2 IFREMER  MODEL  CODE  FROM  CCRS 


%  function  sigmaO  =  if remer (V, phi , theta) ; 

%  calculate  radar  cross-sections  for  cmodS  model  wind  retrieval 


g, 

o 

g, 

o 

g, 

o 


C 


X 


algorithm  according  to  IFREMER  off-line  analysis 
inputs 

V  -  the  wind  speed  in  m/ s 

phi  -  the  relative  radar  look-wind  in  degrees ...  this  angle 
zero  whe  nthe  wind  is  blowing  towards  the  radar 
theta  -  the  incidence  angle  in  degrees 
output 


simgaO  -  radar  cross-section  (linear 


=  [-2.437596  -1.567030 
-0 . 039714  0 .410275 

-0 . 026785  0 . 130832 

0 .192522  -0 . 027294 

0.003595  0.073938 

0.012663  0.044884 

-0 . 007864  -0 . 001635] 

=  (theta  -  36) ./19; 


0 . 370824 
0 . 184862 
-0 .466007 
0 . 032854 
0 . 12845 
0 . 004448 


units ) 


is 


pO  =  1; 
pi  =  x; 

p2  =  (3.*x.^2  -  1)  ./2; 

p3  =  (5.*x.^2  -  3) .*x./2; 

alpha  =  c(l)  +  c(2).*pl  +  c(3).*p2  +  c(4).*p3; 

beta  =  c(5)  +  c(6)  . ’^pl  +  c(7)  .’^p2; 

delta  =  c(8)  +  c(9).’^pl  +  c(10).’^p2; 

bO  =  10.^  (alpha  +  beta . sqrt  ( V  -  delta)); 
bl  =  c(ll)  +  c(12).’'pl  +  c(13).’'V; 

b2  =  c(14)  +  c(15).^pl  +  c(16).^p2  +  (c(17)  +  c(18).^pl  +  c  ( 1 9 )  .  ^p2 )  . 

+  c  (20)  .  ^V.  "^2; 

sigmaO  =  b0.'^(l  +  bl . '*'cos  (phi  .  ^pi/180  )  +  b2  .  ^cos  (2  .  ^phi  .  ^pi/180  )  ) 


C.  INTERNET  SITES: 


CCRS/RADARSAT,  http://www.ccrs. nrcan.gc.ca/ccrs/tekjd/radarsat/rsate.html/ 
ENVISAT/ASAR,  http://envisat.estec.esa.nl/instruments/asar/ 

ALOS/PALSAR,  http://alos.nasda.go.ip/index-e.html 
JPL/NASA,  http://southport.ipl.nasa.gov/desc/imagingradarv3.ltml 
PIRATA,  http://www.brest.ird.fr/pirata/pirataus.html 
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